Entanglement in a Dimerized Antiferromagnetic Heisenberg 

Chain 



Xiang Hao and Shiqun Zhu* 

School of Physical Science and Technology, Suzhou University, 
Suzhou, Jiangsu 215006, People's Republic of China 

Abstract 

The entanglement properties in an antiferromagnetic dimerized Heisenberg spin- 1/2 chain are 
investigated. The entanglement gap, which is the difference between the ground-state energy and 
the minimal energy that any separable state can attain, is calculated to detect the entanglement. It 
is found that the entanglement gap can be increased by varying the alternation parameter. Through 
thermal energy, the witness of the entanglement can determine a characteristic temperature below 
that an entangled state can be obtained. The entanglement detected by the energy can provide a 
lower bound for that determined by the concurrence. If the alternation parameter is smaller than 
a critical value, there is always no inter-dimer entanglement in the chain. 
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I. INTRODUCTION 



The quantum entanglement is considered as key resources of quantum information pro- 
cessing [1, 2]. The property of entanglement plays an essential role in understanding and 
quantifying the physical systems [3-9]. In recent years, some useful measures of entangle- 
ment have been proposed. The entropy of entanglement is used to qualify the entanglement 
of pure states [3] while the entanglement of formation is one measure for mixed states [4]. 
The entanglement in optical systems has been theoretically analyzed and observed [5-9]. 
For further developing the experimental detection, the separability criterion [10, 11] was 
suggested. On this basis, the thermal energy [12, 13], the magnetization [14] and the sus- 
ceptibility [15] have been used as entanglement witness [16] to detect the entanglement. If 
the expected value of the witness is negative, an entangled state could be obtained. These 
measurements can provide an intuitive way to evaluate the quantum entanglement. With 
the help of these measurements, it is possible to find some entangled systems suitable for 
quantum computation and quantum communication [17]. Recently, much interest has been 
focused on the entanglement in solid-state systems [18-27]. The spin chain is one kind of en- 
tangled systems in condensed matters. The entanglement has been observed experimentally 
in Heisenberg spin chains [18]. The dependence of entanglement on the external magnetic 
field [19] and anisotropy [20] has been investigated. The frustration [21, 22] and arbitrary 
spin-s Heisenberg chain [23, 24] were considered. Some interesting effects were discussed 
concerning the relation of entanglement with correlation [25] and quantum phase transition 
[26, 27]. The dimerized spin chain is another essential sort of spin models in many real solids. 
It is often used as a model to explain thermodynamical properties of many substances [28]. 
It is interesting to note that there is the same structure of phase diagrams in the antiferro- 
magnetic dimerized Heisenberg spin- 1/2 and uniform spin-1 chains [29]. The model belongs 
to quasi one-dimensional magnets of spin-ladder systems with even number of coupled spin 
chains [30]. Therefore, it is necessary to evaluate and detect the quantum entanglement in 
alternating spin chains. 

In this paper, the entanglement in an antiferromagnetic dimerized Heisenberg spin chain 
is investigated. In Sec. II, the energy is introduced as a witness for the detection of the dimer 
entanglement. The characteristic temperature for the presence of the dimer entanglement 
is derived. The entanglement gap is introduced to evaluate the entanglement for a chain 
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with large number of spin dimers. The upper bound to the entanglement gap is derived. 
In Sec. Ill, the relation of the entanglement witness and concurrence is deduced. As an 
example, the four-spin dimerized chain with two dimers is analyzed in detail. The limit 
cases are discussed when the alternation parameter equals to zero or one. A brief discussion 
concludes the paper. 

II. WITNESS FOR DIMER ENTANGLEMENT 

The Hamiltonian of the dimerized Heisenberg spin chain can be written as [28, 29] 

L 

H = jj2^-(-ms t -s l+1 (i) 
i=i 

— * -. 

where Si = is the ith spin vector. The number of spins L is even and the periodic 
boundary condition L + 1 = 1 is assumed. At temperature T, the thermal equilibrium state 
is p = exp(—H(3)/Z where Z = Tr[exp(— Hf3)] is the partition function and (3 = 1/kT. For 
the convenience, the Boltzmann constant k is assumed to be one. The values of the exchange 
coefficient J > and J < correspond to the antiferromagnetic and ferromagnetic cases 
respectively. The parameter 5 denotes the alternating ratio of exchange interactions. There 
is an equivalent expression of the Hamiltonian [28, 29] H = j' S2i+aS 2 i- S 2 i+i) with 

J = J'^p and 5 = where a is the alternation parameter. Such dimerized Heisenberg 
spin chain is schematically illustrated in Fig. 1. The spin chain is constructed by the 
dimers of the number d = L/2. Thus, the Hamiltonian can be described as H = Yld=i ^d- 
In Fig. 1, the spin dimer is labelled by an elliptical box. The solid line represents the 
dimer interaction j' and the dash line denotes the inter-dimer interaction aj' . In the 
following discussions, the antiferromagntic case of f — 1 and < a < 1 is considered. 
It is easily found that the ground state at a = is a tensor product of a singlet state. 
When a ^ 0, the ground state cannot be expressed as a tensor product of each dimer 
state. For the simplest case of d — 2, the ground state at a — 1 can be written as l^o) — 
^=[(|1100) + |0011) + |1001) + |0110>) — 2(|1010> + |0101))]. Here |1), |0) are assumed to be 
the eigenstates of the pauli operator a z with the eigenvalues ±1. Therefore, it is clear that 
there exists the transition of the ground state |^> ) if ol is varied from zero. 

On the basis of the separability criterion, there exists a minimal separable energy E sep 
[12, 13]. That is, E sep is the minimal energy that any separable state can arrive at. For 
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the antiferromagnetic dimerized Heisenberg chain, the pure separable state for the minimal 
energy can be analyzed by the standard symmetry methods [24, 31] and expressed as \ip sep ) = 
Ylf =1 (\0)i-\i)i)<S>(\0) i+1 +\l)i +1 ). The minimal separable energy is E sep = (ip se p\H\ip sep ) = 
— L ( 1 + a ) . For the detection of the dimer entanglement, the entanglement gap gE, which is 
the difference between the minimal separable energy E sep and the ground-state energy Eq, 
can be introduced. The upper bound of the entanglement gap can be derived. The lower 
bound to the ground state energy of the Hamiltonian 

H = f y^fla-i • S 2i + aS 2i ■ S 2i+1 ) (2) 

can be estimated. By the variation principle, one can show that the ground state energy of 
the Hamiltonian H can be written as 

E (H) > £ (#i) + a£ (#2). (3) 

— * — * — * — * 

Where E (Hi) is the ground state energy of the Hamiltonian Hi = ^A^i-i-S^) = SiS 2 + 

S3S4 + , while E (H 2 ) is the ground state energy of the Hamiltonian H 2 = J2(S 2i S 2i+ i) = 

S 2 S 3 + S4S5 + Since 

( \ o\ 

0-120 
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(4) 



one has £ (#i) = ^0(^2) = -f-L. Thus for the ground state E (H) > -\L{a + 1). The 
expression of the entanglement gap gE and its upper bound can be given by 

g E = ^ T ^<\(a + l). (5) 

If the entanglement gap is nonzero, the entanglement could be detected experimentally below 
a certain temperature. To some degree, the entanglement gap can be regarded as a useful 
indicator of the existence of entanglement. The entanglement at a finite temperature can 
be detected more easily if the entanglement gap is large [12]. For the dimerized chain with 
a large number of dimers, the ground-state energy E can be numerically calculated. The 
dependence of the entanglement gap on the alternation parameter is plotted in Fig. 2. It 
is shown that the entanglement gap is decreased with the increase of the alternation a and 
reaches a minimal value at about a = 0.7, and then increased again until a = 1.0. Here 



the dimer interaction is chosen as J' = 1. It is found that the entanglement gap qe can be 
increased when a ^ 0.7. Meanwhile, the entanglement witness at finite temperature T can 
be introduced as [12, 13] 

W = E- E sep (6) 

where E = J2d Tr(p^d) i s the energy at the thermal state p. If the value W < 0, an entangled 
state is obtained. The value of W > means that it is a separable (unentangled) state. The 
entanglement witness W as a function of the temperature T and the alternating parameter 
a is plotted in Fig. 3. Figure 3(a) is a three-dimensional plot of the entanglement witness 
W as a function of a and T. From Fig. 3(a), it is seen that the value of W is increased 
to a maximum and then decreased when a is increased. The value of W is increased from 
negative to the positive when the temperature T is increased. That is, the thermal states 
are changed from the entangled states to the separable states as T is increased. The contour 
of W as a function of the characteristic temperature T c and the alternation parameter a 
is illustrated in Fig. 3(b). The value of T c is located on the contour map of W — 0. The 
contour map of W < corresponds to the entangled thermal states. It is found that both 
entanglement-detecting methods are almost equivalent to each other. 



III. RELATION OF ENTANGLEMENT WITNESS AND CONCURRENCE 



Although some of the real solids is composed of a large number of spins, many properties, 
like thermal and magnetic properties can be efficiently studied by the model of small number 
of spins. It is necessary to analyze the entanglement properties of the dimerized chain with 
small number of dimers. In the Hilbert space {|11), |10), |01), |00)}, the reduced density 
matrix of one dimer p d can be written as 

/ u d \ 

w d t d 

t d w d 

\ u d ) 



Pd = 



(7) 



where the elements are expressed as u d = \ + ^f£i[p d cr • a], w d — \ — ^Tr[pd<? • <?] and 
t d = |Tr[p rf <7 • a]. The concurrence of the system is given by C = max{0,2Ai — £V Aj} 
where Aj are the square roots of eigenvalues of the matrix R = p(a y <g> a y )p(a y <g> a y ) in 
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decreasing order and Ai is the maximum one [4]. From the definition of the concurrence 
C, the dimer entanglement C of the alternating Heisenberg chain can be calculated as 
C — 2max{0, \u d \ — \t d \} — \ max{0, — Tr[p d a- a] — 1}. Here Tr(p d a ■ a) < is the correlation 
function which is monotonically increased with the temperature T [23]. According to Eq. 
(6), the witness W can also be expressed by the correlation function 

W = ^r^[l + Tr(p d a.a)} (8) 

Thus, the relation of the entanglement witness W and concurrence C can be expressed as 

AW 

C>max{0,-— - — }. (9) 
(1 + a)L 

Here the equality can be achieved if the alternation parameter a — or 1. The item of 
max{0, — } denotes the entanglement witnessed by the thermal energy. From Eq. (9), 
the detection of the entanglement can offer a lower bound for the concurrence C. 

In the limit case of a = or a = 1, the dimerized Heisenberg spin chain is reduced to 
a single dimer or an isotropic Heisenberg chain. For the case of a — 1, the concurrence 
for any one dimer is just the entanglement between two nearest neighboring spins. It can 
be written by C i>i+1 (a = 1) = \ max{0, -K i>i+1 - 1}, where K i>i+1 = Tr[p^ • a i+1 ] is 
the correlation function for any two nearest neighboring spins. According to Eq. (6), 
the entanglement witness can be expressed by W(a = 1) = E + j. Here the minimal 
separable energy is E sep = — j. For such an isotropic spin chain, the thermal energy is 
E(a = 1) = Y^i Tr(p5i • S^+i) = ^K iti+1 . It is easily seen that the concurrence C iji+ i(a = 
1) = max{0, [24]. When a = 0, it is easily seen that there is no inter-dimer interaction, 

i. e., the correlation function i^2i,2?+i(ct = 0) = 0. The correlation function K 2 i-i^i can still 
exist in the chain. The entanglement witness for a = can also be written by W(a = 0) = 
Y^l 2 Tr(pS 2 i-i ■ S 2 i) + f = f (K 2 i-i,2i + 1)- Similarly, the concurrence for a single dimer is 
C 2 i-i t2 i(a = 0) = max{0, — \{K 2i _ 1>2i + 1)} = max{0, ^^}. Therefore, the equality in Eq. 
(9) is obtained when the alternation parameter is in the limit of a = or a — 1. 

As a simple example, the antiferromagnetic dimerized Heisenberg chain with two dimers 
is analyzed. Because the Hamiltonian satisfies [H, £\ Sf] = 0, the total spin number S is 
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conserved. The eigenvalues E t of the Hamiltonian can be written as 

E 1 = E 2 = E 3 = E 10 = E 13 = ei = (1 + a)/2, (10) 
E A = E 9 = E 14 = e 2 = -(1 + a) /2, 
E 5 = E 7 = E 15 = e 3 = (1 - a)/2, 
E % = E s = E 16 = e 4 = -(1 - a)/2, 
E u = e 5 = -(1 + a)/2 + Va 2 - a + 1, 



#12 = ee = -(1 + a)/2 - Va 2 - a + 1, 
The corresponding eigenstates can also be expressed by 

= |1111); |^ 2 ) = |0000); (11) 
|^3) = ^(|01H) + |1101) + |1011) + |1110)); ^4) = |(|0111> + |1101) - |1011) - |1110)); 

|V> 5 ) = |(|0111> - |1101) + |1011) - |1110)); |Ve) = ^(|0111) - |H01) - |1011) + |1110)); 

\rh) = ^(10011) - 11100)); |^s) = ^(11001) - 10110)); = ^=(|1010> - |0101»; 

|Vio> = 4=(|1010) + 10101)) + |1001) + 10110) + |1100) + 10011)); 
V6 

= "7===i=====M|1100) + 10011)) + y + (|1010) + 10101)) + 11001) + 10110)]; 

v v 1 x + 1 y+) 

IV12) = -==1=== Ml 1100) + 10011)) +y_(|1010) + 10101)) + 11001) + |0110)]; 
y/2(l + xi +y±) 

IV13) = i(|1000> + 10010) + 10100) + 10001)); \4>u) = i(|1000> + |0010) - 10100) - 10001)); 
\fa 5 ) = -(|1000) - 10010) + 10100) - |0001)); \ip 16 ) = -(|1000) - |0010) - |0100) + |0001>) 

where the parameters of x± and y± are given by x± = — [(1 — a) ± \J a 2 — a + 1] and 
y± = —a ± y/a 2 — a + 1. At finite temperature T, the correlation function K = Tr(p^cT • a) 
can be calculated as 

K = -3+^(20e- ei/3 +6e- e2/3 +12e- e3/3 +6e- e4/3 + ~ _ e~ e5l3 + " + 9 e~^) (12) 

Z y l+x 2 _+y 2 _ l + x 2 + + y\ ' v ' 

From Eqs. (8)-(12), the concurrence C and the entanglement witness W of one dimer can 
be easily calculated. 

It is also known that there is inter-dimer interaction in the antiferromagnetic dimerized 
Heisenberg spin- 1/2 chain. The concurrence C 2 i,2i+i between two spins 2i and 2i + l denotes 



the inter-dimer entanglement. It is induced from the the inter-dimer exchange coupling aJ'. 
The reduced density matrix between any two spins has the same structure as that given by 
Eq. (7). Therefore, the entanglement C it j between any two spins i and j is also dependent 
on the correlation function -fQj, which can be expressed by 

C iJ = ±max{0,-K iJ -l}. (13) 

The inter-dimer correlation function is K 2 i t 2i+i = Tr[p 2 i,2j+i^'o1> where P2i,2i+i is the reduced 
density matrix for inter-dimer two spins. It is clearly seen that there is the inter-dimer 
entanglement if the correlation function K 2 i t2 i+i < — 1- For such spin chain, the correlation 
function between any two spins is increased with the increase of the temperature [23]. Thus, 
the inter-dimer entanglement of the concurrence C^j+i is always decreased from that of 
the ground state. When the correlation for the ground state K% ij2i+1 > —1, there is no 
inter-dimer entanglement at any temperature T. As an example, for the case of L — 4, 
the ground state is given by the state | ^12 ) - The inter-dimer correlation function for the 
ground state is obtained by K®^ = —3(1 — 1+x 2 +y 2 )■ In the limit of a — 0, K$ 3 = 0, there 
is no inter-dimer entanglement. When the value of a is increased from zero, the value of 
inter-dimer correlation function is also decreased from zero. If a > 0.5, K% 3 < — 1, the inter- 
dimer entanglement exists. In the limit of a — 1, K® 3 = —2, the inter-dimer entanglement 
reaches the maximum value. The inter-dimer entanglement of the concurrence 6*2,3 can be 
calculated by the correlation function of K$ 3 . By numerical calculations, it is found that 
there is no inter-dimer entanglement when the alternation parameter is below a critical value 
of a c . The critical alternation parameter a c is plotted in Fig. 4(a) as a function of spin 
number L. It is seen that the value of a c increases and then saturates to a constant of 
a s c ~ 0.78 at L = 12. To show the saturation effect, the characteristic temperature T c of 
W = is numerically calculated as a function of the alternation parameter a and is plotted 
in Fig. 4(b) when the number of spins L is varied. It is seen that the value T c of L = 10 
(the bullet) keeps almost the same value as that of L = 12 (the triangle) for relatively large 
value of the alternation parameter a. 

It is very interesting to note that the entanglement between two spins depends on the 
separable distance — in the antiferromagnetic dimerized Heisenberg spin-1/2 chain. This 
can be illustrated by the correlation function K^ - between two spins 1 and j for the ground 
state. For the chain with spin number L = 12, the correlation function K®j = Tr[p° j<7 • a] 
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is numerically calculated and plotted in Fig. 5 as a function of separable distance \j — 1|. 
Here p° j i s t ne reduced density matrix between spins 1 and j for the ground state. In Fig. 
5, the dashed horizontal line denotes the value of K\ j = —1. Below the horizontal line, 
the entanglement exists. From Fig. 5, it is clear that the values of the correlation function 
is K®j > —1 when the separable distance is \j — 1| > 2. According to the relation of the 
correlation function and the concurrence, the concurrence Cij is always zero for \j — 1| > 2. 
In the Heisenberg chain, only K® 2 is smaller than —1 while the others K^^ j>2 ^ > — 1. 
The entanglement always declines from the ground state at T = 0. Therefore, at any 
temperature T, there is no entanglement between two spins i and j when the separable 
distance is \j — i\ > 2 in the antiferromagnetic dimerized Heisenberg spin- 1/2 chain. 



IV. DISCUSSION 



The entanglement in the antiferromagnetic dimerized Heisenberg spin-1/2 chain is ana- 
lyzed. The entanglement gap qe can be increased if the alternation parameter a / 0.7. If 
there is a large number of dimers in the chain, the characteristic temperature T c determined 
by the entanglement witness is decreased to a constant value. The entanglement can be 
detected below a certain temperature in real solids. The relation of the witness and the 
concurrence is also derived. The energy as the entanglement witness can provide a lower 
bound for the dimer entanglement. In the limit of a — or 1, both of the witness and 
the concurrence are equivalent to detect the entanglement. It is also found that there is no 
inter-dimer entanglement when the alternation parameter is smaller than a critical value. 
These results may be helpful for the further study of solid-state quantum communication 
and computation. 
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Figure Captions 

Fig. 1. 

The antiferromagnetic alternating Heisenberg spin-1/2 chain can be described by spin 
dimers. The spin dimer is labelled by an elliptical box. The solid line represents the dimer 
interaction and the dash line denotes the inter-dimer interaction. 

Fig. 2. 

The entanglement gap gE is plotted as a function the alternation parameter. The minimal 
gap is located at about a = 0.7. 
Fig. 3. 

The witness W is plotted as a function of the alternation parameter a and the temperature 

T. 

(a) . The three dimensional plot of W as a function of a and T. 

(b) . The contour plot of W on the plane of (T — a). 
Fig. 4. 

The characteristic alternation parameter a c and temperature T c are plotted. 

(a) . a c is plotted as a function of the number of spins L. 

(b) . T c is plotted as a function of a with L = 6(o), 10(«), 12(A). 
Fig. 5 

The correlation function of the ground state is plotted when the spin separable 
distance \j — 1| is increased from 1 to L — 1. Here spin number L = 12 and the alternation 
parameter a is varied from a = 0.3(o) to a — 0.8(«) to a = 1(A). The dash line denotes 
the value —1. 
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